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We study the Brownian motion of a classical particle in one-dimensional inhomogeneous environ- 
ments where the transition probabilities follow quasiperiodic or aperiodic distributions. Exploiting 
an exact correspondence with the transverse-field Ising model with inhomogeneous couplings we 
obtain many new analytical results for the random walk problem. In the absence of global bias 
the qualitative behavior of the diffusive motion of the particle and the corresponding persistence 
probability strongly depend on the fluctuation properties of the environment. In environments with 
bounded fluctuations the particle shows normal diffusive motion and the diffusion constant is simply 
related to the persistence probability. On the other hand in a medium with unbounded fluctuations 
the diffusion is ultra-slow, the displacement of the particle grows on logarithmic time scales. For the 
borderline situation with marginal fluctuations both the diffusion exponent and the persistence ex- 
ponent are continuously varying functions of the aperiodicity. Extensions of the results to disordered 
media and to higher dimensions are also discussed. 



I. INTRODUCTION 

The Brownian motion is perhaps the best understood 
stochastic process in classical physics, both in homoge- 
neous environments and disordered media. The study 
of the diffusion problem in inhomogeneous environments 
is physically motivated by transport processes (molecu- 
lar diffusion, flow lines in a porous medium, electrical 
conduction) on one hand and the relaxational properties 
of disordered systems (random magnets, spin-glasses) on 
the other (see Refs. [1,2]). 

In the presence of asymmetric transition rates, i.e., 

■ when the the probability per unit time u> r , r ' for a parti- 
cle to jump from site r to site r' is different from w r ',r, 

, the disorder strongly modifies the behavior of the diffu- 
' sive motion in d < 2 dimensions. In one space dimen- 
sion, where the effect of disorder is most pronounced, 

■ the diffusion is ultra-slow and the averaged mean-square 
displacement grows on a logarithmic time scale [3] : 

[(X 2 (t))] av ~\n 4 t. (1.1) 

Another type of inhomogeneity is provided by fractal lat- 
tices, either regular or random, such as percolation clus- 
ters, in which the Brownian motion has been intensively 
studied under the name "the ant in the labyrinth" [4,5,2]. 
Also much work has been devoted to the clarification of 
diffusion processes in the presence of hierarchically dis- 
tributed energy barriers [5,6], a problem which is related 
to relaxation processes in disordered systems [7] . 



In the present work we study the Brownian motion in 
inhomogeneous environments where the transition rates 
are asymmetric and distributed according to quasiperi- 
odic or, more generally, aperiodic rules. As a related 
earlier work, one may mention an investigation of the 
Brownian motion on the two-dimensional Penrose lat- 
tice, where normal diffusive behavior has been found [8]. 
Here we mainly concentrate on one-dimensional aperiodic 
systems. Besides its mathematical interest, the present 
study is also physically motivated since artificial multi- 
layer systems with controlled distributions of the atomic 
layers may now be grown by molecular beam epitaxy [9] . 
When particle transport has different time scales for the 
motion parallel and perpendicular to the layers, respec- 
tively, a one-dimensional diffusion process perpendicular 
to the layers is, in principle, observable. 

The study of cooperative phenomena in quasiperiodic 
and aperiodic systems is an intensive field of research. 
One may mention phase transitions and critical phenom- 
ena in Ising and other magnetic models, percolation, self- 
avoiding walks, etc. Aperiodic structures, which interpo- 
late between periodic and random systems, may or may 
not influence qualitatively the properties of a cooperative 
process. Concerning the critical behavior of aperiodic 
magnetic systems, a relevance-irrelevance criterion has 
been proposed [10,11], which is an extension of the well- 
known Harris criterion for disordered systems [12]. The 
vast amount of exact results about the critical properties 
of aperiodic quantum Ising chains [13-17] and related 
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aperiodically layered two-dimensional Ising models [18] 
are all in accordance with this criterion. Apcriodicity 
may also change to second order a transition which is of 
the first order in the pure system, as was demonstrated 
recently in a numerical study [19]. 

In the present work we show that a relevance- 
irrelevance criterion, similar to that of magnetic systems, 
can be formulated for the Brownian motion, which is then 
checked against exact results obtained for different onc- 
dimcnsional aperiodic environments. 

The paper is organized as follows. In Section II we in- 
troduce the basic notations and quantities (drift velocity, 
diffusion constant, persistence probability) and present 
the relevance-irrelevance criterion for one-dimensional 
aperiodic environments. In Section III an exact cor- 
respondence between the random walk (RW) and the 
transverse-field Ising model (TIM) is presented in one 
dimension, which is then used in Section IV to obtain an- 
alytical results for irrelevant, relevant and marginal ape- 
riodic environments. Our results are extended to higher 
dimensions in Section V and discussed in the final sec- 
tion. 



II. FORMALISM AND 
RELEVANCE-IRRELEVANCE CRITERION 

We consider a one-dimensional RW with nearest neigh- 
bour hopping, characterized by transition probabilities 
Wij±i for a jump from site i to site i ± 1. The time evo- 
lution of Pi(t), the probability for the particle to be on 
site i at time t, is governed by the Master equation: 



1 T1 N 
1 ~ rin=l e n 



Vd 



1 lln=l tn ^ N Z^n=l Z^i=l fc n+l t n+2 • ■ ■ t n+i 



Defining the control parameter as 



1 N 

JvE lne - 



(2.3) 



(2.4) 



the drift velocity is zero for 5 = 0, whereas for a small 
bias, such that SN <C 1 , the drift velocity is proportional 
to 5: 



1 



N N 



\f2 ^2 ^2 e n+l 6 n+2 ■ ■ ■ e n+i ■ ( 2 ' 5 ) 



N 2 



n—1 i—1 



One can similarly calculate the diffusion constant D, 
which in the zero bias case, 5 = 0, is simply given by: 



D(5 = 0) = D . 



(2.6) 



Before going to analyze the diffusive behavior of aperi- 
odic walks we consider another quantity of interest, the 
persistence probability P per (L,t), which is the probabil- 
ity that the walker starting at site i = 1 does not cross 
its starting position until time t. Here the length scale 
L in the definition is set by the presence of an adsorbing 
site at i = L + 1, thus wl+i,l = 0. Due to this ad- 
sorbing site the persistence probability has a finite long 
time limit, lirn^oo P per (L, t) = p per (L), which can be 
expressed as [21]: 



dPi, 
dt 



(2.1) 



The transition probabilities are generally nonsymmctric, 
here we suppose that their ratio is given by: 



Wj,i+1 



= eR fi 



(2.2) 



where R > is the amplitude of the inhomogeneity 
(R = 1 corresponds to the homogeneous lattice) and the 
integers fi are taken from an aperiodic or quasipcriodic 
sequence. For the sake of simplicity in the following we 
take u>i+i,i = w^- = 1. 

The aperiodic chain may be replaced by a periodic ap- 
proximant of period N, such that e, = e i+ M [20,10]. The 
aperiodic system is recovered in the limit N — > oo. Exact 
expressions for the drift velocity Vd and the diffusion con- 
stant D obtained by Dcrrida for an arbitrary distribution 
of the transition rates [20] can then be used to treat the 
aperiodic case. 

With our notations, the drift velocity of the aperiodic 
model may be written as: 



Pper 



(2.7) 



i=l j=l 



It is easy to see that, in the thermodynamic limit, 
liniL^oo Pper {L) = p per plays the role of an order pa- 
rameter: It is nonvanishing for S > only. For the ho- 
mogeneous system with R = 1 in (2.2): 



horn 



(2.8) 



whereas at the critical point in a finite homogeneous sys- 
tem: 



p h p °™(L,d = 0) = 



1 



L + l 



(2.9) 



P P er(t) = limL^oo P per (L, t) is the usual persistence 
probability for an infinite system. 

Next we analyze the expressions of the basic quanti- 
ties in Eqs. (2.5) and (2.7). The qualitative behavior for 
any periodic system, i.e., when N is finite, is the same in 
the vicinity of the critical point: The diffusion constant 
Dq in Eq. (2.5) is finite and the persistence probability 
in (2.7) has a linear 8- or L~ ^dependence. 
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For aperiodic systems, which are obtained in the limit 
N —* oo, one has to consider the products of the tran- 
sition rates V n (L) = Ylf =1 e n+i = exp[A„(L)], where 

A n (L) = X)f=i me "+i measures the fluctuations in the 
"energy landscape" when In a is related to the heigth of 
the i-th energy barrier in an activated diffusion process. 
Here we differentiate between three possibilities: 

i) For bounded fluctuations in the transition rates, i.e., 
with a flat energy landscape, the products V n (L) stay 
finite at the critical point, thus the diffusion constant D 
remains finite and the linear ^-dependence of the persis- 
tence probability is maintained. Thus this type of aperi- 
odicity is irrelevant and the aperiodic Brownian motion 
keeps the same critical properties as the homogeneous 
one. 

ii) For unbounded fluctuations in the transition rates, 
when the energy landscape is rough, so that A„ (L) <~ L n 
with a wandering exponent Q > 0, the products V n (L) 
which appear in Eqs. (2.5) and (2.7) are divergent at 
the critical point when L — > oo. Consequently, the dif- 
fusion constant D vanishes and the persistence proba- 
bility displays a nonlinear 5-dcpendcncc. We conclude 
that an aperiodic environment with unbounded fluctua- 
tions is a relevant perturbation for the Brownian motion. 
As we show later, diffusion in a relevant aperiodic envi- 
ronment is ultraslow, i.e., the mean-square displacement 
grows on a logarithmic time-scale, as for the Sinai model 
in Eq. (1.1). 

iii) In the borderline case, the wandering exponent of 
the environment is = 0, thus fluctuations in the energy 
landscape grow logarithmically. In this marginal situa- 
tion the products of the transition rates V n {L) have a 
power low L-dependence and the different physical quan- 
tities also display power law singularities. For this type 
of anomalous diffusion the drift velocity depends on the 
bias as: 

v d ~ S T , \8\ « 1 , (2.10) 

with r > 1, whereas the mean square displacement grows 
like 

(X 2 (t))~t^, (2.11) 

with ip < 1. Finally, the persistence probability is in 
general characterized by algebraic singularities: 

Pper ~5 x , P per (t)~t- e . (2.12) 

The critical exponents defined above are not indepen- 
dent. Using results of the next section in Eqs. (3.7) 
and (3.9), one can show that the relevant time scale t 
and the bias 6 are related by t^ 1 <~ v^S ~ S 1+T . For the 
length scale, one has L <~ v^t <~ S^ 1 . These relations lead 
to the scaling laws 

r=\~^ X=|e, (2.13) 

which evidently hold for normal diffusion with r = ip = 
X = 1 and 6 = 1/2. Thus two critical exponents are 



enough to describe the behavior of the anomalous diffu- 
sion. 

To obtain the characteristic quantities for relevant and 
marginal environments one has to solve the Master equa- 
tion in (2.1), which amounts to solve the eigenvalue prob- 
lem for the transition matrix or Focker-Planck (FP) op- 
erator, 

Kv k = \ k v k , ulK^ul^k, (2.14) 

where the matrix elements take the form = Wij 

for i ^ j and (.1/); , = — J2j w i,i- Ah 1 the physical prop- 
erties of the model can be expressed in terms of the left 
and right eigenvectors u k and y_ k , respectively, and the 
eigenvalues X k , which are nonpositivc. For example, the 
probability Pij(t) that the walker starting at t = at 
site i arrives on site j at time t is given by: 

p iM = ^WfeWwfe(j)exp(A fe i) . (2.15) 
k 

The relevant time-scale of the problem is related to the 
inverse of the largest nonzero eigenvalues and the dynam- 
ical properties of the RW are connected to the scaling be- 
havior of the eigenvalues of the FP operator at the top of 
the spectrum. Considering a large finite system of size L, 
under a change of the length scale by a factor of b > 1, 
such that L' = L/b, the eigenvalues at the top of the 
spectrum with fc<L are expected to transform as: 

A' fc = 6»»A fc , (2.16) 

where y\ is the scaling dimension of the gap. From 
Eq. (2.16), the finite-size behavior of the eigenvalues 
Afe(L) ~ L~ Vx follows, thus the time and length scales are 
related through t ~ L yx . On the other hand from (2.11) 
L 2 ~ t^, thus the diffusion exponent is given by: 

V = -. (2.17) 

y\ 

For relevant aperiodic environments, the leading eigen- 
values have a stretched exponential finite-size depen- 
dence, thus the diffusion exponent is formally zero. 

In the Master equation formalism, the persistence 
probability P per (L,t) can be calculated by putting ad- 
sorbing sites at i — and i = L + 1. Then P per (L,t) = 
Pi,L+i(t), which in the large-i-limit is just the first com- 
ponent of the zero-mode left eigenvector p per (L) = ui(l), 
as given by Eq. (2.7). 



III. CORRESPONDENCE BETWEEN THE 
RANDOM WALK AND THE ISING MODEL IN A 
TRANSVERSE FIELD 

First we rewrite the eigenvalue problem of the FP op- 
erator in (2.14) in terms of the components of the right 
eigenvector v k (i) as: 
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Wi-i,i v k (i - 1) - (wis-i + Wi, i+1 ) v k {i) + Wi+i.i v k (i + 1) = A fe v k (i) , (3.1) 
and consider a finite system of size L, i.e., we put w .i = Wl+i,l = 0. Then we introduce the new variables: 

TT^±I] , (3.2) 

in terms of which the eigenvalue problem is transformed into: 

(wi_ M w M _i) 1/2 v k (i - 1) - {wis-i + w hi+1 )v k (i) + (w i+M w M+ i) 1/2 v k (i + 1) = XkV k {i) , (3.3) 



which corresponds to a real symmetric eigenvalue prob- 
lem Y,j T i]^k{j) = Awfc(i) with Tij = Tji. Consequently 
the eigenvalues Afe of the FP operator are real. 

Next we show that the eigenvalue problem in (3.3) for- 
mally appears in the free-fermion representation of the 
transverse-field Ising spin chain, with couplings J, and 
transverse fields hi, described by the Hamiltonian: 

L-\ L 

i=l i=l 



Comparing Eqs. (3.6) and (3.3), one can notice that they 
can be cast into the same form with the correspondences: 

T . / a/2 

hi (wi, i+ i) 1/2 ^ 7 \ 

-A fe . 

Thus there is a mathematical equivalence between the 
RW in an inhomogeneous environment and the TIM with 
the corresponding inhomogeneous couplings, as described 
in (3.7). One can however go further and show that there 
are several physical quantities which are closely related in 
the two problems. 

First let us consider the persistence probability 
Pp er (L), which is calculated with adsorbing boundary 
conditions ti>o,i = u>l+i,l = as p per (L) = ui(l)v\(L + 
1) = ui(l)a^ 1 a L+ iVi(L + 1) = [vi(L + l)] 2 . Now using 
the correspondences in (3.7), the equivalent TIM with 
ho = Jl = has a fixed spin at i = 0, whereas the other 
end of the chain at i — L+l is free [23]. The surface mag- 
netization of the chain measured on the i = L + 1-th spin 
is given by m s (L) = $i(L + 1) [23,24], thus using (3.7) 
we have a relation between the surface magnetization of 
the TIM and the persistence probability of the RW: 

[m s (L)} 2 ^ p per {L). (3.8) 



I 

where af and erf, are Pauli matrices at site i. The Hamil- 
tonian H can be transformed into a free-fermion model 
by standard techniques [22]: 

ff = £A fc (r,+ % -l/2). (3.5) 

k 

Here the (r/k) are fermion creation (annihilation) op- 
erators and the excitation energy Afe is the solution of 
the eigenvalue equation: 



(3.6) 
I 

We note that this relation has already been mentioned in 
Rcf. [21]. 

Let us now analyze the relation between the energy 
scales in the two problems, Afe and Afe, respectively. 
First one may notice that, according to the last equa- 
tion in (3.7), the eigenvalues Afe of the FP operator are 
nonpositive, as they should on physical grounds. 

For the periodic TIM with Ji = h and hi = hi+N = 

1 /2 

he i ' , the energy of low- lying modes, close to the critical 
point 5 — 0, are given by a perturbation calculation [10] 
as: 

Alc^v^ + ql), (3.9) 

where q k ~ 1/N denotes the wavenumbcr corresponding 
to the smallest excitation energy. Comparing the ex- 
pression of the sound velocity v s in Ref. [10] to D in 
Eq. (2.5), we obtain another correspondence: 

v 2 s <^=> D . (3.10) 

Taking the infinite periodic approximant limit, TV — > oo, 
we obtain the following important result: if in an envi- 
ronment the diffusion is normal, i.e., the diffusion con- 
stant Do is finite, then the phase transition of the cor- 
responding TIM with inhomogeneous couplings is in the 



I 

J 4 _i/i,_i$fe(« - 1) + (J 2 _! + h 2 t )<S> k {i) + Jihi$ k {i + 1) = k 2 k $ k {i) ■ 
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Onsager universality class. On the other hand, if the dif- 
fusion is anomalous, i.e., the diffusion exponent is ip < 1, 
then the phase transition of the TIM is not of Onsager 
type. In the marginal case the phase transition in the 
TIM is characterized by an anisotropy exponent (or dy- 
namical exponent) z > 1, which describes the finite size 
scaling behavior of the gap A& <~ L~ z for k <C L. Ac- 
cording to Eq. (3.7), z is related to the diffusion exponent 
through: 




On the other hand, in the relevant case, A^ has a 
stretched exponential L-dependence, thus the dynamical 
exponent is formally infinite. 

Finally, we present a useful estimate for the smallest 
(in absolute value) nonzero eigenvalue of the FP op- 
erator, X m im by transforming a related result for the 
TIM [16,23]. Having a large finite chain of length L with 
reflecting b.c. at i = 1 (u>i,o = u>o,i = 0) and adsorb- 
ing b.c. at i = L + 1 (w l +i,l = 0, ^ 0), then 
the leading finite-size behavior of X min (L) is connected 
to the product of the persistence probabilities at the two 
ends of the chain as: 

L 

X min (L) ~ - Pper (L)p per (L) He; 1 , (3.12) 

i=l 

where the persistence probability at site L is given by: 

P P er(L)= (l + ErK-; J ' ( 3 - 13 ) 

whereas p per (L) is given by Eq. (2.7), however for L — 1 
sites. 



where poc is the asymptotic density of the perturbed 
transition probabilities. Then the fluctuation of the "en- 
ergy landscape" in a large system of size L is character- 
ized by: 

L 

A(L) = J^lne, = n L \nR + Llne c 

»=i 

= \nR{n L - L Poo ) ~ In R L n . (4.2) 

Here fi is the wandering exponent of the aperiodic se- 
quence [25] which is easily obtained for aperiodic substi- 
tutional sequences. 

Working with a finite alphabet A, B,..., such that, via 
subsitution, A — > S(A), B — > S(B), etc, the sequence fi 
in (2.2) is obtained by starting with one of the letters, 
iterating the substitution process and finally replacing 
the letters by digits (or groups of digits) and 1. For a 
two-letter sequence, one can directly proceed with subsi- 
tutions on and 1. 

The fluctuation properties of a sequence can be de- 
duced from its substitution matrix with entries riy giv- 
ing the numbers of letters i — A,B,... in S(j) (j — 
A,B,...) [25]. The wandering exponent fi involves the 
two largest eigenvalues in modulus of this matrix, ^ and 
fi2, and reads fl = ln |/i2|/ hi/Ui. 

We now consider different environments associated 
with specific aperiodic sequences. Environments with 
bounded (ft < 0), unbounded (fi > 0) and marginal 
(O = 0) fluctuations are treated separately. Most of the 
results can be obtained by adapting the analytical meth- 
ods developed for the aperiodic TIM , using the corre- 
spondences given in (3.7), (3.8), (3.10) and (3.11). 



A. Aperiodic environments with bounded 
fluctuations 



IV. DIFFUSION IN APERIODIC 
ENVIRONMENTS 

Here we consider different one-dimensional aperiodic 
environments and study the corresponding diffusive be- 
havior by analytical methods. As we mentioned before, 
the relevance-irrelevance of the perturbation is connected 
to the fluctuation properties of the aperiodic environ- 
ment. 

To be more specific, we consider the behavior at crit- 
icality (5 = 0) of the logarithm of the transition rates 
lne^, which is related to the heigth of the i-th energy 
barrier, Ae i: in an activated diffusion process. When the 
global bias vanishes, e takes its critical value e c such that, 
according to (2.4), 

L 

e c = R~ p °° , poo = lim n L /L , n L = V fi , (4.1) 

L^oo ^— ' 
i=l 



1. Quasiperiodic (Fibonacci) environment 

Quasipcriodic lattices can be generated in several dif- 
ferent ways, among others by the well known cut-and- 
project method. Here we use the following algebraic def- 
inition for a one-dimensional quasiperiodic sequence: 



i 













where [x] is the integer part of x and lu > 1 is irrational. 

The Fibonacci sequence can be generated by the sub- 
stitution rule -> 010, 1 -» 01 [26] starting with 0. 
When read from left to right, it corresponds to (4.3) with 
lu = (\/5 + l)/2, the golden mean. 

The sequence in (4.3) and the corresponding quasiperi- 
odic environment, as defined in Eq. (2.2), have bounded 
fluctuations since the wandering exponent in (4.2) is 
fi = —1. Consequently, the diffusion constant D in (2.5) 
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is finite. As shown in Appendix A, it can be expressed 
in closed form using the methods of Ref. [10] as: 

/ lni? \ 2 , A , 

It is interesting to note that D does not depend on the 
value of the irrational parameter u>. 

The persistence probability in (2.7) can be also evalu- 
ated analytically using the techniques of Ref. [27] (see Ap- 
pendix A). For small bias there is a linear J-dependcnce: 

P P er(S) = Rl ~ 1/UJ -j^[ S > 0<S«1, (4.5) 

where the prefactor depends on u>. 

For the other end of the system, except for the first 
digit which is irrelevant, the sequence, read from right to 
left, is given by: 



i + lj 




i + u — 1 


[ _ 




LO 2 



(4.6) 



which leads to the persistence probability: 

p per (<5)=ir 1+1 /- T ^ T <5, 0<<5«1. (4.7) 

One may notice a simple relation between the diffusion 
constant and the persistence probabilities, 

P P er(S)p per (S)=D Q S 2 , (4.8) 

which is valid for any value of the irrational parameter 

U). 



2. Thue-Morae sequence 

The binary Thue-Morse sequence [28] is generated 
by the substitutions — > 01 and 1 — > 10, leading to 
0110100110010110... The wandering exponent is ft = 
— oo thus the sequence has bounded fluctuations. Conse- 
quently, the diffusion constant is finite and, as shown in 
Appendix A, can be calculated along the lines of [10] as: 

^^ U/ 4+ Vv 4 ) 4 - (4 - 9) 

The persistence probability is linear for small bias and 
given by: 

in agreement with the result of Ref. [13], when properly 
translated. 

Since the Thue-Morse sequence is reflection symmet- 
ric, the persistence probability is the same at both ends of 



the system, p per (S) — P per (^) so that the relation in (4.8) 
is again fulfilled. 

Analyzing the results obtained for environments with 
bounded fluctuations we are led to the following conclu- 
sions: 

i) Predictions of the relevance-irrelevance criterion are 
fully satisfied: The diffusion constant is finite and the 
critical exponents (r = ip = \ = 1) take the same values 
as in homogeneous environments. 

ii) In both examples, the diffusion constant D is in- 
variant under the transformation R — > 1/R. This is in 
agreement with the fact that Dq(R) is maximal for the 
homogeneous system, i.e., at R = 1. 

iii) For symmetric sequences, the persistence probabil- 
ities p per and p per are equal, whereas, for asymmetric 
sequences, they satisfy the relation p per (R) = p per (l/R). 



B. Aperiodic environments with unbounded 
fluctuations 



The aperiodic environments display unbounded fluc- 
tuations when the wandering exponent > 0. Conse- 
quently, the diffusion constant in (2.5) vanishes in the 
infinite periodic approximant limit and the persistence 
probability in (2.7), like the surface magnetization in the 
TIM, has an anomalous behavior. 

To obtain a qualitative picture for the behavior of the 
RW in such an environment, we estimate the leading 
eigenvalue of the FP operator A m i„ in a finite system of 
size L, taking into account that the relevant time-scale 
of the process is given by t <~ A m J n . 

According to Eq. (3.12), the size dependence of X m in 
for a sequence with unbounded fluctuations is primarily 
determined by the product n^T] 1 e i~ 1 ~ exp(— constL n ). 
Thus, at criticality, the leading eigenvalue has a stretched 
exponential size-dependence: 

Xmin ~ exp(-const L n ) . (4.11) 

This implies an ultra-slow diffusion process, which takes 
place on a logarithmic time-scale: 

(X 2 (t)) - (lnt) 2 /°. (4.12) 

In a random environment, with fl — 1/2, Eq. (4.12) cor- 
responds to the Sinai diffusion in (1.1). 

Next, we consider the persistence probability p per and 
analyze its properties for an exactly solvable environment 
with unbounded fluctuations. 
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1. Rudin-Shapiro environment - persistence probability 

The Rudin-Shapiro sequence [28] is obtained via sub- 
stitutions on pairs of digits, 00 -> 0001, 01 -> 0010, 
10 — ► 1101 and 11 — > 1110, such that starting on 00, one 
generates the sequence 0001001000011101 . . . The wan- 
dering exponent of the sequence, £1 = 1/2, is the same as 
for a random environment. 

In the critical situation, the persistence probability can 
be exactly calculated using the methods of [14] and the 
correspondence in (3.8). One obtains completely differ- 
ent behaviors for R > 1 and R < 1: 

i) For R > 1, the persistence probability remains finite 
for an infinite system at criticality: 

R 1 ' 2 

Pper = 1 ~ 1 - W 2 + R ' (4 - 13) 

thus the transition in the persistence properties is dis- 
continuous as the bias S changes sign. 

ii) For R < 1, the persistence probability has an 
anomalous size dependence at the critical point: 

Pper ~ exp(— const L 1 / 2 ) , (4-14) 

which goes to zero in the infinite-size limit. 

Reading the sequence from the other end, amounts to 
exchange the digits 1 <-> 0, which is equivalent to the 
transformation R — > 1/R. Thus the persistence proba- 
bility of the critical Rudin-Shapiro environment has very 
different properties for the two ends of the chain: While 
it stays finite at one end, it vanishes at the other. At this 
point, one may ask how the persistence probability be- 
haves if the system starts at an arbitrary position i along 
the sequence. One can answer by using the correspon- 
dences in Section III and the known results for the sur- 
face magnetization of the Rudin-Shapiro TIM model [29]. 
Typically, i.e., with probability one, the critical point 
persistence vanishes, as in (4.14). However, there is a 
fraction of starting points, the so-called "rare events", 
Prare ~ L~ ea ", for systems of size L, where the persis- 
tence probability is finite p per = 0(1), as in (4.13). Then 
the average persistence, which is obtained by averaging 
over the starting positions, scales as the probability of 
the rare events, thus 

[p P er(L)] av ~ Prare(L) ~ L~ @ ^ , (4.15) 

where [. . .] ot , denotes an average over the starting posi- 
tion. The same scenario remains valid for other aperiodic 
sequences with unbounded fluctuations. Translating ex- 
act results for the TIM [29] into the RW language, the 
persistence exponent <d a v is connected to the wandering 
exponent Q, of the sequence via: 

e„„ = i-n. (4.16) 

To summarize, the diffusion process in an aperiodic en- 
vironment with unbounded fluctuations is anomalous: i) 



The average displacement grows logarithmically in time 
and ii) The persistence probability is not self-averaging, 
its averaged dependence involves the fluctuation expo- 
nent n of the environment. 



C. Aperiodic environments with marginal 
fluctuations 

Aperiodic environments with marginal fluctuations are 
characterized by a wandering exponent SI = 0. The crit- 
ical behavior of the TIM with marginally aperiodic cou- 
plings is nonuniversal and several critical exponents vary 
continuously with the strength of the aperiodic pertur- 
bation. From the correspondences presented in Section 
III, it follows that the diffusive behavior of the RW in 
this type of environment is also anomalous. Generally, 
the diffusion constant D vanishes and both the diffusion 
exponent, ip = ip(R) < 1, and the persistence exponent, 
X = x(R) > 1 are continuous functions of the asymmetry 
parameter R. 

For aperiodic environments with marginal fluctuations, 
the scaling behavior of the eigenvalues of the FP opera- 
tor at the top of the spectrum can be obtained exactly by 
a renormalization group (RG) transformation, as intro- 
duced for the TIM in Ref. [15], later applied to different 
sequences in [16] and generalized in [17]. 

The essence of the method is a decimation procedure 
such that, after one step, a fraction 1/6 of the original 
lattice sites are left. The master equation for the renor- 
malized system can be cast into the same form as the orig- 
inal one, provided a finite set of appropriate new parame- 
ters are introduced into the original equation (parameter 
space). The linearized transformation at the fixed point 
corresponding to A* = gives the gap exponent y\ as 
indicated in Eq. (2.16). To construct explicitly the RG 
equations, we refer to the related work on the TIM in 
Refs. [15-17]. The same type of RG procedure can be 
used to calculate the persistence exponent, but one can 
also simply deduce it from a finite-size-scaling analysis 
of Eq. (2.7) at criticality as shown in Appendix B for a 
specific sequence. 

In the following we present results translated from 
Refs. [15,16] for several environments with marginal fluc- 
tuations. 



1. Period- doubling environment 

Using the substitutions — > 01 and 1 — > 00 and 
starting with 0, one generates the period-doubling se- 
quence [30] 0100010101000100 . . . which, apart from the 
last digit, is symmetric and has a vanishing wandering 
exponent. 

The diffusion exponent in the period-doubling environ- 
ment is given by: 
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In 2 



In (JJV6 + R-l/6) ' 



(4.17) 



whereas the persistence exponents, which are the same 
at both ends due to symmetry, read: 



In (i?V6 + jj-i/6) 
X = X = — 



In 2 



(4.18) 



as shown in Appendix B. 



To summarize, the diffusion in marginally aperiodic en- 
vironments is anomalous, both the diffusion exponent tp 
and the persistence exponents \i X are continuous func- 
tions of the parameter R, they however satisfy the scaling 
relation 



x(R) + x(R) 



which is a consequence of (3.12). 



(4.24) 



2. Paper-folding environment 

The paper-folding sequence can be generated by a re- 
current folding of a sheet of paper [28]. It corresponds 
to the two-letter substitutions 00 -> 0010, 11 -> 0111, 
10 -> 0110 and 01 -> 0011. Starting with 00 one arrives 
at 0010011000110110 . . .. This environment has the same 
(R ^ type of asymmetry as for the Rudin-Shapiro 

sequence, if one forgets the last, irrelevant, digit. The 
wandering exponent vanishes. 

The diffusion exponent in the paper-folding environ- 
ment is 



In 2 



In (flV4 + #-1/4) 



(4.19) 



The persistence exponents are different at the two ends 
of the system and are given by: 

lnfl + i?- 1 /*) ln(l + i?V2) 
X= V — n X= y -• (4-20) 



In 2 



In 2 



3. Hierarchical environment 

Here we consider the Huberman-Kerszberg hierarchi- 
cal environment [31], where the positions i of the digits 
fi satisfy the relation: 



i = 2 f '(2l + l), 1 = 0,1,2 



1 ' ' ' 5 



(4.21) 



thus the sequence starts as 0102010301020104 . . .. The 
diffusion problem in the same environment with symmet- 
ric transition rates has already been studied before [6] . 
In the nonsymmetric case the diffusion exponent is 



In 2 



In + fl-i/2) 



(4.22) 



The persistence exponents are different at the two ends 
of the system and are given by: 



In (1 + R~ 
In 2 



X 



ln(l + fl) 
In 2 



(4.23) 



One can also easily solve the diffusion problem for gen- 
eralized hierarchical environments following the solution 
of the corresponding TIM in Ref. [16]. 



V. RELEVANCE-IRRELEVANCE CRITERION IN 
HIGHER DIMENSIONS 

The relevance-irrelevance criterion of Section II can be 
generalized for a d-dimcnsional environment, where the 
nonsymmetric transition rates w r y are perfectly corre- 
lated in (d — D) dimensions. Thus they vary in D < d 
dimensions and the fluctuations are characterized by a 
wandering exponent f2. 

In our considerations the basic role is played by the 
displacement probability P(L), which measures the frac- 
tion of walks that have moved to a distance L from 
their starting position during time t ~ L 2 , which is 
the characteristic time-scale for a homogeneous medium. 
In a one-dimensional homogeneous environment with a 
weak uniform bias < S u = ln(w^/w^) <C L , the 
displacement probability in the unfavorable direction is 
P U (L) ~ exp(— const 5 U L), thus exponentially small. 

In the absence of a global bias, an inhomogeneous (ran- 
dom or aperiodic) environment does not generally favor a 
net displacement of the particle. But due to the fluctua- 
tions in the transition rates, the motion can be effectively 
biased locally. We now estimate the effective force or av- 
erage local bias 6(L) inside a domain Ut, of linear size L. 
First, generalizing the expression in (4.2), we calculate 
the accumulated value of \n(w T y /w r > , r ) in the domain 
as: 



ML) = Yl H 

reu L 



w r 



>/Wr>, T ) - lni? L° n . 



(5.1) 



Then the averaged local bias along the inhomogeneous 
directions is given by: 



5{L ) ~ ^ ~ L -D ( i-n) 



(5.2) 



Thus the displacement probability, in analogy with the 
uniform case, is given by: 

P{L) - cxp[-const 8(L)L] - exp [-const L 1 --^ 1 -^) . 

(5.3) 

Now, depending on the sign of the exponent 

<t> = 1 - D(l - Q) , (5.4) 
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the behavior of the displacement probability and, conse- 
quently, the diffusion properties are different. 

i) For <fr < 0, i.e., for CI < 1 — 1/D, the displacement 
probability has no exponential size dependence, thus one 
has the same type of diffusive behavior as for homoge- 
neous systems with zero bias. Consequently, the envi- 
ronment does not modify the normal diffusive motion of 
the particle and therefore this type of perturbation is ir- 
relevant. 

ii) On the other hand, for <fi > 0, i.e., for £1 > 1 — 1/D, 
the displacement probability decays exponentially with 
L^, thus this type of environment is relevant for the 
diffusive properties. The relation between the relevant 
time-scale t, which is proportional to the characteristic 
number of steps needed to have a displacement L, and 
the length scale L is obtained as: 

t ~ P(L)- 1 ~ cxp (const L - Dil - n)+1 ^ . (5.5) 

Thus the mean square displacement, which is propor- 
tional to L 2 , grows on a logarithmic scale as: 

(X 2 (t)) ~ (lni)i-°( 2 i-») . (5.6) 

iii) The borderline case, <p — 0, corresponds to the 
marginal situation where the L-dcpcndcnce in the expo- 
nential of Eq. (5.3) can be logarithmic, leading to nonuni- 
versal diffusive behavior, as observed in Section IV for 
one-dimensional systems. 

We note that the above relevance-irrelevance criterion 
is in complete agreement with the exact results we ob- 
tained above for d = D = 1. For example Eqs. (4.12) 
and (5.6) coincide in this case. On the other hand, in the 
case of uniform disorder with d = D and Q = 1/2, the 
borderline dimension predicted by Eq. (5.4) is d* = 2, 
which is in agreement with the results of RG investiga- 
tions [32]. 

VI. DISCUSSION 

In this paper we have studied the scaling properties 
of the Brownian motion in inhomogeneous environments 
where the transition rates are asymmetric and their vari- 
ation follows some quasiperiodic, aperiodic or hierarchi- 
cal rules. It has been shown that the diffusive motion of 
a particle in such type of environments can be anoma- 
lous and a relevance-irrelevance criterion has been for- 
mulated, which allows to predict the different scenar- 
ios. In one dimension, we have obtained many exact re- 
sults which all are in agreement with the abovementioned 
relevance-irrelevance criterion. In these one-dimensional 
calculations, we have made use of an exact mathematical 
correspondence between the Master equation of the RW 
and the eigenvalue problem for the energy of the free- 
fcrmionic excitations of an inhomogeneous TIM. This 
correspondence has been exploited to obtain relations be- 
tween different physical quantities in the two problems. 



The analytical results previously obtained for the TIM 
have been translated into exact results for the diffusion 
problem. 

At this point we note that there is another model of 
statistical physics, the directed walk, which is also closely 
connected to the TIM. As was shown in Ref. [15] the 
scaling properties of the directed walk are connected to 
the eigenvalues of Eq. (3.6) at the top of the spectrum. 
Thus the three problems (Ising model, random walk and 
directed walk) are inherently related, the complete solu- 
tion of any of those contains the necessary information 
about the properties of the two others. In particular, one 
single RG transformation describes the scaling properties 
of the three models: the fixed point at A = for the TIM 
governs the critical properties of the Ising model and that 
of the Brownian motion, whereas the fixed point at the 
top of the spectrum is connected to the properties of the 
directed walk. 

Next one can show that the parametrization of the 
transition probabilities, = 1 below Eq. (2.2), does 

not affect the conclusions of the paper. The relevance- 
irrelevance criterion in Section II is evidently unaffected 
by this restriction, and in the marginal situation the 
nonuniversal exponents %, \ and ip are also insensitive 
to this parametrization. For the persistence exponents it 
follows from the fact that in Eqs. (2.7) and (3.13), only 
the ratio of the transition rates appears. From Eq. (4.24), 
or more generally from the relation in Eq. (3.12), the 
same conclusion is reached for the exponent tp. This re- 
sult can be also obtained by analyzing the structure of 
the systematic RG-technique of Ref. [17]. 

The properties of the Brownian motion in a relevant 
aperiodic environment are in many respect similar to 
those in a disordered media. In the critical situation, 
i.e., in the absence of global bias, the diffusion is ul- 
traslow in both cases and the corresponding relations 
in (1.1) and (4.12) are analogous. In the off-critical situ- 
ation, however, there is an important difference between 
the diffusive properties in the two environments. In a 
disordered media for small enough global bias, such as 
< S < S + , the diffusive motion of the particle is anoma- 
lous and the average displacement of the walker grows as 
a power law, [{X(t)}} av ~ t^, with a (^-dependent ex- 
ponent < yu(<5) < 1 [33]. This regime is in complete 
correspondence with the Griffiths-McCoy phase of the 
random transverse- field Ising spin chain [21]. As shown 
very recently in Ref. [34] , this anomalous diffusion regime 
is absent for relevant aperiodic environments. 
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APPENDIX A: DIFFUSION CONSTANT AND 
PERSISTENCE PROBABILITY FOR THE 
FIBONACCI AND THUE-MORSE 
ENVIRONMENTS 

Using (4.1), the expression for the diffusion constant 
D in (2.5) can be rewritten at criticality as: 

N N 

— - — W f- n - 1 R~ 

Do ~ N 2 2^2^ e c * 
n—1 i—1 

N N 

= ^££^- 9 " +i > ( A1 ) 



where g i = n i - ip^. 

Similarly, the persistence probability in Eq. (2.7) is 
given at criticality by: 



1 



Pper(L) 



i + £)ir*. 



(A2) 



»=i 



For the aperiodic environments, in the limit N — > oo 
and L » 1, respectively, the powers of R in Eqs. (Al) 
and (A2) can be replaced by their averaged values. 

i) For the Fibonacci sequence, gi + l/u) is the fractional 
part of (i+l)/u>, which is uniformly distributed over [0, 1] 
for uo irrational. Thus, in the limit iV — > oo, the average 
in (Al) can be replaced by an integral, 



1 

A) 



R 9 dg [ R~ 9 dg, 
I) Jo 



(A3) 



which leads to the result given in Eq. (4.4). 

In the same way, for the persitence probability, 
Eq. (A2) can be rewritten as: 



1 f 1 
= LR 1 ^ / R- 

Pper J 



g dg, 



(A4) 



which is evaluated in (4.5) using the correspondence be- 
tween S and L~ x . 



ii) For the Thue-Morse sequence, the average can be 
simply performed by noticing that fn-x+fn = 1, 01 and 
10 appearing with the same probability , thus gn = 
with probability 1/2 and (fe+i = ±1/2, each with prob- 
ability 1/4. Then considering the different parity combi- 
nations for n and i in (Al), one obtains the expression 
given in (4.9). In the same way, (A2) leads to (4.10) for 
the persistence probability. 



APPENDIX B: FINITE-SIZE SCALING 
CALCULATION OF THE PERSISTENCE 
EXPONENT FOR THE PERIOD-DOUBLING 
ENVIRONMENT 

For the period-doubling sequence, fzi = 1 — fi and 
./Wi = 0, so that: 



ri2i = n 2 i+i = y^j2k=i-n l . 
fe=i 

On a system with size AL, one has: 



(Bl) 



4L 



(B2) 



and, splitting the sum into odd and even parts: 
S , 4i(e, R) = 

2L 2L-1 
= 1 + e- 1 + (T 2i Rr n ^ + e -2i-lfl-n2 i+ i 



i=l 
2L 



i=l 



2L-1 



= 1 + e- 1 + Y(t 2R )~ lRn * + ^ J2 (t 2 R)~ l R nt 
i=l i=l 

~(l + e- 1 )5 2L (e 2 i?,i 1 '- 1 ). (B3) 
Iterating one step further, the last relation leads to: 

S 4L (e, R) = (l + e" 1 )(l + e- 2 R- r ) S L (e 4 R, R) . (B4) 
At the critical point, e c = i? -1 / 3 , this gives: 

S4 L (e c ,R) = (R 1/6 + R- 1/6 ) 2 S L (e C7 R). (B5) 

Since p per ~ S x ~ L~ x , one immediately recovers the 
persistence exponent given in Eq. (4.18). 
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